Damping in 2D and 3D dilute Bose gases 
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Damping in 2D and 3D dilute gases is investigated using both the hydrodynamical approach 
and the Hartree-Fock-Bogoliubov (HFB) approximation . We found that the both methods are 
good for the Beliaev damping at zero temperature and Landau damping at very low temperature, 
however, at high temperature, the hydrodynamical approach overestimates the Landau damping 
and the HFB gives a better approximation. This result shows that the comparison of the theoretical 
calculation using the hydrodynamical approach and the experimental data for high temperature 
done by Vincent Liu (PRL 21 4056 (1997)) is not proper. For two-dimensional systems, we show 
that the Beliaev damping rate is proportional to and the Landau damping rate is proportional 
to for low temperature and to T for high temperature. We also show that in two dimensions 
the hydrodynamical approach gives the same result for zero temperature and for low temperature 
as HFB, but overestimates the Landau damping for high temperature. 



INTRODUCTION 



The experimental realization of Bose-Einstein conden- 
sation (BEC) in magnetically trapped akali atoms 
provides a good tool to study the properties of 3D dilute 
Bose gases. Furthermore, with the anisotropic of new 
traps [1, , one can further confine condensate atoms 
in a quasi-two-dimensional regime § . Most theoreti- 
cal work (for 3D and 2D systems, see Review 0, 0: re- 
spectively) has focused on the dynamics of condensates 
and the zero-temperature behavior, which can be ob- 
tained by solving the non- linear Gross-Pitaevskii (GP) 
equation. However, the finitc-temperature behavior has 
still remained difhcult to study, where experimcnts have 
shown damping of the condensates modes in the presence 
of a significant noncondensate component [l^ . Damp- 
ing mechanism associated with collective excitations of 
Bose condensed atoms interacting with a non-condensed, 
thermal component is not well understood and still rep- 
resents a challenging problem in theoretical physics. The 
damping of collective modes can have various origins. 
There are two distinct contributions to the total decay 
rate 7 = 7_b + 7l: One arises at T = O from the pro- 
cess of decay of a quantum of excitation into two or more 
excitations with lower energy. This mechanism was first 
studied by Baliaev |Tl|in 3D uniform Bose superfluids 
and is known as Balieav damping 75. At finite tempera- 
ture, a different mechanism of damping (known as Lan- 
dau damping, 7^) comes from the process of one quantum 
of excitation decays due to coupling with transitions as- 
sociated with other elementary excitations and occurring 
at the same frequency. Landau damping is not associated 
with thermalization process and can be well described in 
the framework of mean field theory T^, [11, [13. The 
subject of Landau damping in dilute BECs has been ex- 
plored by several authors. Landau damping in an uni- 
form Bose gas at low temteratures was first investigated 
by Popov [15] Hohcnberg and Martin [l^, while at high 
temperatures was first investigated by Sz'epfalusy and 
Kondor The relevance of Landau damping to ex- 



plain experimental data of trapped Bose gas was pro- 
posed by Liu and Schieve jl^ and developed by Liu [l^ 
using the Popov's hydrodynamical approach . On the 
othe hand, Pitaevskii and Stringari [l^l investigated Lan- 
dau damping in a weakly interacting uniform as well as 
non-uniform Bose gas by means of semi-classical theory. 
They showed that for the uniform Bose gas, it reproduces 
known results for both the low temperature asymptotic 
behaviour of the phonon coupling. However, for the high 
temperature, Liu showed higher Landau damping rate 
than those obtained by Szepfalusy and Kondor, while 
the hign-temperature behavior could be reproduced by 
Pitaevskii and Stringari. 

However all investigations of the damping rate has 
been done for a 3D Bose gas. 2D Bose gases are in- 
teresting as their low temperature physics is governed 
by strong long-range fluctuations. These fluctuations in- 
hibit the formation of truc long-range order, which is a 
key concept of phase transition theory in 3D. Thus a 
2D uniform interacting Bose gas does not undergo Bose- 
Einstein condensation at finite temperatures. However, 
this system turns superfluid below the BKT (Berezinski, 
Kosterlitz and Thouless) temperature Tkt [13, ISI . The 
experiment indication of the BKT transition in weakly 
interacting Bose system has even been shown in Ref. [g] . 
Damping in a 2D Bose gas is an open question which was 
recently addressed by the authors for a uniform Bose gas 
[23*1 using the hydrodynamical theory of Popov [15]. In 
this work, we show that the hydrodynamical approach 
actually overestimates the damping rate at high tem- 
peratures, both for 3D and 2D systems and we calcu- 
late the Balieav and Landau damping rates for a 2D 
uniform Bose gas using the semi-classical Hartree-Fock- 
Bogoliubov (HFB) approach. In the limit of low temper- 
atures, the results of this approach is in good agreement 
with that found from hydrodynamical approach. Con- 
trary to earlier work [1^, we show that for the 3D case 
in the high temperature limit, the hydrodynamical ap- 
proach cannot be used to explain the experimental data. 

This paper is organized as foUows. In Sec. |TT]we dis- 
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cuss the relation between the atom-atom interaction and 
the scattering length for 2D and 3D dilute gases. In 
Sec. |lll]we first introduce the hydrodynamical approach 
developed by Popov^l5'|, and then calculate the BeUaev 
damping and Landau damping for 3D and 2D gases. The 
mistake using this approach for high temperature is also 
discussed. In Sec. IIVI the HFB approximation is devel- 
oped to calculate 3D and 2D Beliaev and Landau damp- 
ing. 



II. ATOM-ATOM INTERACTION AND 
SCATTERING LENGTH 



The Standard Hamiltonian of an interacting Bose gas 



is 



d'*rd'^r'V'^(r)V'^(r')?7(r - r')i^{r)i^{r'), 



(1) 



where U is the atom-atom interaction and V^xt 
external potential. For uniform Bose gas, V^xt = 0. The 
true interaction between atoms is very complicated where 
one has to consider the fine structure of atoms. However, 
the scattering process can offer a effective potential to 
simplify the ineraction. In order to do that, one has to 
introduce Green functions for bosonic systems with con- 
densate. The difhculty of doing so arises from the fact 
that the terms containing the odd number of annihila- 
tion operators do not vanish for a Bose gas after averag- 
ing the ground state due to the existence of condensate, 
which unfortunately destroys the hope to apply the nor- 
mal technique of Feynman diagrams to the system. This 
difficulty was successfully resolved by Beliaev pH . [23 |. 
He separated the operators with zero momentum, which 
semi-classically can be regarded as a c-number, and the 
other operators with nonzero momenta. In this way, the 
Feynman diagrams can be used for the Bose gas. 

Beliaev considered a three-dimensional system with 
short-range, central interaction potential with radius 
1/03 and then calculated the renormalized atom-atom in- 
teraction in the presence of the condensate between two 
particles with non-zero momenta, which one should sum 
over all ladder diagrams. In this way, one can obtain 
the renormalized interaction in terms of the s-wave scat- 
tering am plit ude according to the elementary scattering 
theory [tI. l25l. [2^ . Therefore the effective potential can 
be written as 



f/(r) 



(2) 



with the atom mass m, and the momentum dependence 
of the scattering amplitude can be ignored in the low 
temperature limit. In the rcst of the paper we define the 
atom-atom interaction strength 173 as 



53 



m 



(3) 



For two dimensions, Schick foUowed the methods de- 
veloped by Beliaev and examined a two-dimensional sys- 
tems of hard-core bosons with a diameter 02 at low den- 
sity and zero temperature. Unlike the three-dimensional 
systems, where ladder diagrams are independcnt of the 
dimensionless parameter najj, hence, it is natural to take 
it as the small perturbation terms to expand the quan- 
tities. For two-dimensional systems, contributions from 
the ladder diagrams depends logarithmically on na| , the 
dimensionless parameter for 2D systems, but not directly 
on itself. In particular, the renormalized interaction 
is proportional to 1/ln 1/7102: 



47r 



92 



min (l/na2 



(4) 



Schick concluded that the l/lnl/na2 plays a role of a 
small parameter in the two-dimensional dilute systems, 
and other quantities, like damping rate in this paper, can 
be expanded in terms of it. 



III. HYDRODYNAMICAL APPROACH 

In the low temperature and low energy limit, Popov 
[lil i developed a hydrodynamical approach to find an 
effective Hamiltonian for a nonideal Bose gas. In or- 
der to do that, one has to separate the order parameter 
over rapidly and slowly oscillating field, and the hydrody- 
namical Hamiltonian can be obtained by integrating the 
functional over rapidly oscillating field. The theory de- 
scribes then the hydrodynamical Hamiltonian in terms of 
two slowly varying fields : phase (t>{x) and density fluc- 
tuation 7t{x) — n{x) — no with no- the density of the 
ground state. Here the four-dimension Euclidean space 
X = (x, r) is used with the imaginary time r. The hydro- 
dynamical action for a d— dimensional nonideal Bose gas, 
according to Popov, can be written in the form (notice 
that h — 1 throughout the paper) 



rl3 

dr 
'o 

1 d^p 
'29^ 



d'^x 



i— — —TTOr(p 



1 dp 



dfidn 
1 d'^p 



(V0)= 



2m d/i 
8mno 2m 



(5) 

with the atom mass m, the pressure of a homogeneous 
system p, the chemical potential /x and the atom density 
71. The fields (j) and tt are periodic in imaginary time 
T with the period (3 ~ l/{kBT). For very low tempera- 
ture, as long as the non-condensate part can be neglected 
compared to the condensate, the pressure 71) at zero 
temperature can be a good approximation. Therefore we 
can use the expression of a weakly interacting dilute gas 
a,s p = fin — ^71q, where gd is the atom-atom interaction 
related to the scattering length. It foUows that 



d^p 
d^dn 



1; 



dp 



d^P 



0; 



d'^p 



-9d, 



(6) 
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<(p(kj <p(-k)> «P(kj TT(-k}> <w(k) 0(-k)> 



<TT(k) TT(-k)> 




vertex 



FIG. 1: Green's functions and vertex 




and the action 



dTd'^x i iird. 



takes the form 



(np + tt) 
2™ 



9d 2 (V7r)^\ 



(7) 

The action contains aU quadratic functions except the 
term l/2m7r(V(/))^, considered as an interacting poten- 
tial. The Hamiltonian can be derived from the effective 
action ([7]) as 



,a , w 2 
d X I -^n-v 



9d, 



where the field of velocities is defined as v = l/niWcj). 
This Hamihonian is consistent with one particular real- 
ization of the Landau hydrodynamical Hamiltonian. 

Fourier transforming the ficlds 4> and tt, the effective 
action IH]) can be writtcn in the form 



1 



2 ^ 



d'^k—k^(j){k)(l){-k) + 2uj^(l){k)n{~k) 



{dd 
1 



Amn, 



-)7r(fc)7r(-fc) 



o 



(9) 



3V ^ 



2m 



(/.(fci)0(fc2)7r(fc3) 



where k is the vector (k, ioj^) and the Matsubara frequen- 
cies coi, = 2-Kv/ j3 with integers v. From the action of the 
Fourier transformation ^ one can extract the important 
information needed for the perturbation calculations us- 
ing diagrammatic technique. First of all, the free Green's 
functions is defined as foUows 



^o^''^ - i (7r(fc)0(-fc))o {n{k)n{-k))o, 



(10) 



where (• • • )o denotes the expectation value of fields cal- 
culated only with the quadratic action. From the action 
©, the inverse of the free Green's function can be found 
as 

"o Jj2 

m 



(11) 



Therefore, 



Go(fc) = 



ffd+k /(imno) 



4mno 



(Tio/m)k^ 
w2+£^(k) 



(12) 



FIG. 2: One-loop diagrams for the self energy. 



where 



6(k) 



k2 

2m 



c2k2 



(13) 



with c = gdno/m. We represent the relation be- 
tween the free Green's functions and the Feyman di- 
agrams in Fig[T] The cubic term of the action 
^ is known as phonon-phonon interaction in the 
low tcmpcraturc region, giving rise to a vertex of 

(5''(ki + k2 + k3)(5,,j+^2+^3,o^^^7;r^ represented by the 
last diagram of Fig. [1] 

The exact Green's function has to involvc the phonon- 
phonon interaction, given by the Dyson cquation G{k) = 
Go{k) + Go(fc)E(fc)G'(fc), where S(fc) represents the self- 
energy matrix. The low-frequency spectrum of coUective 
modes can be obtained by the poles of the exact Green's 
function as 

detG-\k) = det[Go\k)) - E(/fc)] = O, (14) 

through the analytical continuation iuji, = lu + irj (r; = 
O"*") after the Matsubara frequency sum. The complex 
frequency lu — E — ?7(k) represents the energy spectrum 
E and the damping rate 7. Neglecting the matrix E, 
the zero order approximation for the Eqs. (fTi|) gives the 
square of the Bogoliubov energy spectrum: 



E^ 



Zm 



.(k)^ 



(15) 



When the phonon-phonon interaction is considered, the 
imaginary part appcars in the spectrum. Fig. [5] shows 
the onc-loop diagrams for the self energy S. The contri- 
bution to the imaginary part of the spectrum is given by 
the last five Figs. [2jb) -[2l(f). The damping rate can be 
obtained from Eqs. ^ and ([M]) [ioll as 



7(k) 



i9d 
nok^ 



Amno 
ImE^7r(k, UJ + irj) 



)ImE00(k, o; -f iri) 



(16) 



ReE07r(k, UJ + ir]). 
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Replacing ([15]) into Eqs. (fT6|) and calculating the dia- 
grams, we obtain 7(k) ~ 7B(k) +7i(k), where 7^ is the 
Beliaev damping as 



(k-k')^(k-k0^e(k')£(k) 
2TOnok2k'^e(k - k') 



7B(k) = ^,^^,,1 I d'^fc'<5(e(k) - 6(k') - e(k - k')) 

[/°(e(k'))-/°(-6(k-k'))] 

^ (k-k')(k-(k~kO)e(k) 
2mnok^ 
k^(k'-(k-k'))2e(k')£(k-kO 
4TOnok'^(k-k')2e(k) 
^ (k'.(k-kO)(k.k06(kO | 

and 7l is known as Landau damping as 



(17) 



= 2^+^ / d'fc'-^lflk) + e(k') - 6(k + k')) 



[/"(6(k'))~/"(6(k + k'))] 

[ <k) 

I 2mno 



k'^(k- (k + k'))2e(k + k') 



+ 



k2(k + k')2e(k') 
(k + k')2(k-k')2e(k')l , (k-k')(k- (k + k'))e(k) 



k2k'2e(k + k') 



mrtok^ 



k2(k' ■ (k + kO)2e(k')e(k + k') (k' ■ (k + kO) 
2mnok'2(k + k')2e(k) ^^^^^^^ 
(k-k')e(k') (k-(k + k'))e(k + k') 



(k + k')2 



where the bosonic distribution function 

l/[exp(/3e)-l]. 



(18) 



A. Quantum regime cfc S> fcsT 

At r = O, the Landau damping disappears and the 
Beliaev damping contributes to the damping rate. In the 
Beliaev damping mechanism the momenta of the three 
excitations are comparable, |k| ~ |k'| ~ |k — k'|. Then 
the Eqs. ([nl) yields 



7s(k) 



9c 



d'^fc'lkllk'llk-k'l 



<5(6(k)-e(k')-6(k-k')) 



(19) 



In three dimensions the damping rate for small k {k ^ 
mc) is 



7r^(k) = 



6407rmno 



(20) 



known as Beliaev's resultlllll 



For two-dimensional systems, the Eqs. (IT9)) can be 
written as 



71=' (k) = 2 



dk' 



Ik'|(|k-k' 



sm( 



(21) 



where 6 is the angle between k and k'. The factor two 
in front of the bracket comes from the faet that there are 
two angles corresponding to the energy conservation for 
the Beliaev damping (e(k) — e(k') — e(k — k')): sin0 ~ 

zfc ^2!^*^ ^ ^^]' ''^^'^ Behaev damping rate for a 2-D 
Bose gas has the form 



327rn 



o 



(22) 



This result corrects the wrong result previously given by 
Chung and Bhattacherjee [23| and the factor two will ap- 
pear naturally in the two-dimensional Landau damping. 



B. Thermal Regime cfc > fcsT 

For finite temperature and small momenta such that 
cq <C kgT and cq <^ nogd, the Beliaev damping is much 
smaller than the Landau damping. In three dimensions, 
the damping rate to the lowest order in k can be obtained 
from Eq. (fT8|l as 



e[k) IdirninoksT 



1,5 



(23) 



where 



1 r°° 7 



2z2 



(24) 



'2 2(z2 + l) (z2 + l)2 (z2 + l)2 



with fco — y/rrmog^ and r = ksT/nogs ■ This result was 
first obtained by V. Liu For /csT < mc2^qlM|is 
reduced to the Hohenberg and Martin's result 



Sn^kiksTY 
40mnoc^ 



(25) 



This low temperature limit gives the same result as that 
using the HFB approach, which will be introduced in 
the next section. For high temperature fcsT ^ mc^, 
T{t) ^ 38.735r, and the damping rate is approxiniated 
by 



' ~ 9.648-2 ^ 

e{k) c 



(26) 



with the three-dimensional scattering length 03 — 
mg^/ 4:77. Unfortunately this result is different than that 
investigated by Szepfalusy and Kondor [l7j , which reads 



7l '(fc) _ SnkBTos 
e{k) ~ 8 c ' 



(27) 



Therefore the hydrodynamical approach is no longer cor- 
rect for the high temperature. The reason is that in the 
hydrodynamical Hamiltonian only the slow osciUating 
fields are considered by integrating out the fast osciUat- 
ing fields. For high temperature the fast osciUating fields 
should also be considered to reduce the damping rate. 
We can conclude that the hydrodynamical approach is 
very good for low temperature, however, for high tem- 
perature, other method should be introduced. We will 
discuss that in the next section. 

In Fig. [3] the three-dimensional Landau damping per 
unit energy using the hydrodynamical approach (dashed 
line) and HFB (solid line) is plotted as a function of t. 
Also shown are the asymptotic behavior at high tem- 
perature (the dashed-dot line) and the low temperature 
limit (the dashed-dot-dot line). We can see that the hy- 
drodynamical approach gives very good agreement with 
the low temperature limit for r < 0.5, however, it goes 
too large at high temperature and it does not approach 
the asymptotic value given by Szepfalusy and Kondor. 
Therefore the conclusion made by Liu in Ref. [l^ that 
the results of the hydrodynamical approach can fit the 
experimental data is not proper. 

In two dimensions, the damping rate reads 



where 



sech2(^) 



(28) 



dzz^ 



3 2z2 



(29) 

2 



'2 2(z2 + l) (z2 + l)2 (22 + l)2J- 

In the low temperature limit: ksT ^ mc^, 12(7) - 
-\/67r2T2, therefore the damping coefhcient is given by 



e{k) 8 noc2 • 

In this low temperature regime, the damping rate is pro- 
portional to T^. As far as we know, this quadratic de- 
pendence of the temperature for the damping rate in the 
low temperature is found for the first time in this paper. 

For the high temperature, as the three-dimensional 
case, the hydrodynamic Hamiltonian overestimates the 
damping. In the next section, we will use the 
Hartree-Fock-Bogoliubov approximation to obtain the 
two-dimensional damping at high temperature. 

Fig. ([4]) shows the two-dimensional Landau damp- 
ing rate per unit energy using both hydrodynamical 
approach (dashed line) and HFB method (solid line). 
In this figure the low temperature limit (dashed-dot- 
dot line) and the asymptotic value at high temperature 
(dashed-dot line) are also shown. The hydrodynamical 
result is in agreement with the low temperature limit 
for r < 0.2, however, it will not approach the asymp- 
totic value for high temperature similar to the three- 
dimensional case. 
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FIG. 3: Landau damping rate per unit energy versus r in 
three dimensions. The unit of 7_t/e(k) is ^Jl^rTo. Solid 
black line represents the result obtained by HFB method, 
dashed blue line by the hydrodynamical approach. The high- 
temperature asymptotic behavior and low-temperature limit 
are also shown by dashed-dot green line and dashed-dot-dot 
red line, respectively. 



IV. HARTREE-FOCK-BOGOLIUBOV 
APPROACH 

In this section we represent a semi-classical method: 
Hartree-Folk-Bogolubov (HFB) . We will see that in the 
low-temperature regime this approach is in a good agree- 
ment with the hydrodynamical approach, while for the 
high temperature, on the contrary to the hydrodynamic 
approach, HFB gives a better approximation to the decay 
rate. 

We start with the method by Giorgini fis"] . The grand- 
canonical Hamiltonian of a system with a nonuniform 
external field V^xt('^) reads 



K =H - ^iN = j d''r^T(i.^ t)HoiP{r, t) 

j drV'^(r,t)V'^(r,t)V'(r,i)V'(r,t), 



(31) 



where 



H. 



2m 



+ ^ext(i") - M 



(32) 



and 'ip^T^t) and '4>{Y^t) are the creation and annihilation 
field operators. Since the system is in the regime where 
the condensate exists, we define a time-dependent con- 
densate wave function <I'(r,t) [16] 



$(r,t) = (^(r,t)) 



(33) 



with the average (• • • ) using the grand-canonical Hamil- 
tonian ([?T|) . We have to notice that the Eq. ([55)1 can 
always be used if the condensate exists, however, for a 
homogeneous two-dimensional system, i. e. VQy^i{v) = O, 
the condensate does not exist at finite temperature. In 
this case the long-range order disappears, it remains the 
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quasi-long-range order for a two-dimensional homoge- 
neous Bose gas. That means, though a macroscopic oc- 
cupation number of a single state does not exist, there 
exists a small value of kc in thc momentum space where 
a macroscopic occupation number of the states k < kc 
still forms a quasi-condensate. Therefore the bracket 
in Eq. ([55)1 should count all the states in the quasi- 
condensate. We can see that $(r, t) aUows us to describe 
the oscillating condensate away from the equihbrium. To 
avoid the confusion for the notation, we define here the 
stationary value of the condensate in its equihbrium as 
$o(r) 



$o(r) = (^(r))o, 



(34) 



where (■ • • )o denotes the time-independent average of the 
condensate in its equiUbrium. The particle field can be 
decomposed into a condensate and a noncondensate com- 
ponent 



'ip{r,t) = $(r,i) + Vi(r,t). 



(35) 



By the definition of the condensate ([33]) , the nonconden- 
sate component has to satisfied the condition: 



(36) 



By applying the decomposition ([34| to the grand- 
canonical Hamiltonian, it can be separated to a quadratic 
and a quartic term: K = K2 + ^4, where 



K2= j d-^r (v;^(r,t)ffoV;(r,t) 

+ 2gd|$(r,t) |2 v:;t(r,i)^(r,i) 



-|^$'^^^(r,i)^(r,i)) 



(37) 



and 



Ki^YJ ^)^(^' ^)^(^' (38) 

We are interested in the regime where the condensate 
shghtly differs from the equiUbrium state, that is, 



$(r,t) = $o(r)+<5$(r,t) 



(39) 



with a small fluctuation (5$(r, i). Expanding K2 up to 
the linear term in (5$(r,t), we can rewrite it as K2 — 

Kif'' + K2^\ where K!f ^ is the zero order term of the 
condensate 

-no(r) (^t (r, i)^t (j., t) + ^{v, t)i,{v, t) 



9d^ 

2 



(40) 



with the condensate density no(r) = |$o(r)P, while i^2^^ 
is the linear term in the fluctuation: 



K. 



(1) 



d'^r2grf$o(r)[<5$o(r,t)+<5$S(r, t)] 



4>\r, t)i>{r, t) + .gd$o(r) (5$o(r, i)^^ (r, t)i>\r, t) 



+(5$S(r,i)^(r,t)V'(r,i) . 



(41) 



As for thc quartic term K4, the mean- field decompo- 
sition is first used 



where 



n{r,t) = (V;'f(r,t)'0(r,i)) 
m{r,t) = (^(r,t)^(r,i)). 



(43) 



as the normal and abnormal time-dependent density. Un- 
der the linearization in the fluctuation (|39p . the normal 
and abnormal density are also displaced with a small fluc- 
tuation as 



h{r, t) 
m{r, t) 



n°(r) + Sn{r,t), 
fn^{r) + 5rh{r, t), 



(44) 



expanding around their stationary values n" (r) — 
{ip'<{r,t)ip{r,t))o and m°{r) = {ijj{r,t)^jj{r,t))o. In the 
literature, often referred to so-called Popov approxima- 
tion 



m" (r) = O 



(45) 



has been used in the mean- field treatment. Atually, this 
approximation was never suggested by Popov, as indi- 
cated by Yukalov ^2T\ , but was first proposed by Shohno 
[28j and we will refer it to Shohno approximation or 
Shohno Ansatz in the remaining paper. The Shohno ap- 
proximation is necessary in the treatment because with- 
out that the elementary excitation would have a gap, 
which disobeys the gapless spectrum of the Goldstone 
modes caused by the continuous gauge symmetry break- 
ing in the ground state. However, the using of Shohno 
approximation is still under debate. Several attempts 
have been done to go beyond the Shohno approximation 
(for example, see Ref. [23, H^, H^l). The Popov approx- 
imation is needed in the mean-field factorization (I42p . 
By avoiding factorization, for example, using the per- 
turbation or random-phase approximation to calculate 
the quartic terms, the gapless dispersion can be obtained 
even without Popov approximation. 

Inserting Eqs.(l39l), ^^j, (gH) and the Shohno Ansatz 
(HSI) into the quartic term K4 can be expanded up 



to the first order terms K4 
Sn and S m as 



K. 



(0) 



Ki^^ in fluctuations 



(46) 
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is the zero order term, which represents the coupling to 
the condensate from the quartic term, while the first or- 
der term reads 

^(1) ^9d I ^^d^ (4a(r,i)^t(r,t)^(r,i) 

+ (5m(r,<)'0^(r,t)Vi'f(r, t) + 6m*{r,t)-ip{r,t)-ip{r,t) 

(47) 

Unlike represents the couphng between the fluctu- 
ations of the condensate and the noncondensate parti- 
cles, K^^^ is related to the couphng of the noncondensate 
particles to the normal and abnormal densities. If the 
density of the noncondensate particles is much smaller 
than the density of the condensate, K2^^ is more im- 
portant than k[^\ therefore K^'^ can be neglected and 

In the case of large occupation number of particles in 
the condensate, is much smaller than K2^^ and Kj['\ 
we can use K +Kj[''' as basis to develop a perturbation 
expansion in terms of K2^\ To diagonalize K2^'' + Kj['\ 
one can apply a Bogoliubov transformation 

4,ir,t) = ^w,(r)a,(t)+z;*(r)a](i) 

i 

4'\r,t) = J2u^{r)a]{t)+v,{r)a,{t), (48) 



In order to obtain the decay rate, we have to find 
the time evolution of the fluctuation of the condensate: 
6^{r,t). The equation of motion: 



(53) 



leads to the result 



i— V'(r,i) = HoiPir^t) + gdijHr,tMr,mir,t). (54) 



Inscrting the decomposition (|35p into the equation of mo- 
tion (l54l). it reads 



(55) 



z^$(r,t) = (i/o + .9d|$(r,i)n $(r,0 

+2gd^{r, t)n{r, t) + gd**(r, t)m{r, t) 



Here we assume that the cubic product of the noncon- 
densate contributes very little to the dynamics of the 
condensate, therefore the average value is set equal to 
zero: {i/j^ijjtl;) — 0. The wavefunction in the equilibrium 
can bc obtained by setting d^a/dt = O, which leads to 
the stationary equation 

{Ho + gd [no(r) + 2nO(r)] ) $o(r) - 0. (56) 



By inserting Ea. (j39p and the stationary equation (j56p 
into the equation of motion for the condensate (|55p . the 
equation of motion for the small amplitude 5^ reads 



with the quasi-particle creation and annihilation oper- 
ators aj, a j obeying the Bose commutation relations 
[ai,aj] = 5ij, which gives the normalization condition 
for the functions Uj{r,t), Vj{r,t) as 



J d'^rK(r)u,(r)-<(r)t;,(r)]=5, 



(49) 



Therefore the operator K2^'^ + ifj'^'' can be diagonalized 
if the Bogoliubov-de Genes equations are satisfied: 

Cuj{r) + gnoir)v-iir) = ejUj(r), 

Cvj{r) + gnoir)uj{r) = -ejVj{r), (50) 

where a Hermitian operator is introduced as 

C = Ho + 2gMr) (51) 

with the total density n (r) defined as the sum of the con- 
densate density and normal density in the equilibrium: 
n (r) — no{r) + nP{r). As a result, the grand-canonical 
Hamiltonian (13111 becomes 



K = K2 +Ki = ^ eja\a + K. 



(1) 
2 



(52) 



with the eigenvalues tj obtained from the Bogoliubov-de 
Gennes equations (|50p . 



i^mr, t) = (Ho + 2gdn{r)) 5$(r, t) + .9rf7io(r),5$*(r, t) 
ot 



+25d$o^^(i-, t) + gd$o(r)^m(r, t)- 



(57) 



Applying the Bogoliubov transformation (j48| to Eq. 
([57)) . the Eq. ^ gives the final form: 



z|-<5$(r, t) = {Ho + 25<jn(r)) 5$(r, t) + gdno{r)d<^>* (r, t) 
ot 

+ gd*o(r) ^ {2[u*uj + v*Vj + v*uj]f,j{t) 

+ [2v^Uj + u,uj]g,j{t) + [2v*u* + (t)} , 

(58) 

where fij{t) = (aj(t)aj(i)) and gij = {ai{t)aj{t)) are 
normal and anomalous quasiparticle distribution func- 
tions. 

To calculate the normal and anomalous quasiparticle 
distribution functions using the perturbation Hamilto- 
nian (|52p . we have to use the equation of motion : 



d 



(59) 
(60) 



8 



To the first order, the Fourier transform of and gij at 
the frequency uj is given by 



X $0 [S^i{r,uj)(u,u* + v^v* + v,u*) 



(61) 



d'^r 



LU - {e, + ej)+iO+ 
+ <5$2(r, u;){u*v* + v* u* + v* v*)] , 



(62) 



where S^i{r,oj) and (5<i>2(r, (jj)are the Fourier transform 
of (5$(r,i) and (5$*(r,i): 



(63) 
(64) 



and is the bosonic distribution function 



f ^ : 



(65) 



with f3 — l/ksT. Fourier transforming the equation of 
motion (|58|) and replacing Eqs. ((6T|) and ((62)) into it, we 
obtain the perturbed eigenfrequency: 



and Aij, Bij and Bij are defined as 



A, 



B, 



d'^r^o [5<i>'i{u,u* + v,v* + v^u*) 

J d^r<fo [S<^",{u*v* + v*u* + u*u*) 
+ 6<i>",{u*v* + v*u* + v*v*)] 

d'^r^o [5$? (u., v j + ViUj + UiUj) 



(69) 



The real part of the right-hand sidc (Eq. ([M)) ) gives the 
eigenenergy of the system and the imaginary part teUs 
us the damping coefRcient 7. Using the relation 



X + «0+ X 



iTTS{x), 



(70) 



we can divide the damping rate into two different types: 
one comes from the process that one phonon with the 
frequency loo is absorbed by a thermal excitation 
jumping to another thermal excitation with the energy 
e j — ei+LL}Q. This mechanism is so-called Landau damp- 
ing given by the second term on the right-hand side of 
Eq. (IMI), 

7L = 47r.g^^ I A,,f(/° - /°)<5(c.o + e. - e,)- (71) 

This process happens mostly at finite temperature, it is 
therefore a thermal process. Another kind of decay arises 
from the process of a long wave-length phonon decaying 
into two phonons, as indicated by Beliaev, and it can 
be obtained by the imaginary part of the first term in 
brackets on the right-hand side of Eq. ([66]) : 



(72) 



UJ —ujo + 



45. -/, )-— 7— 



^0 + (ei - Ej) + iO+ 



where the unpeturbed eigenfrequency uJo is obtained from 
the unperturbed RPA equation [3l| 



-g^no -C l\S'fl)-''"[ 5$5 



with the normalization condition 



(67) 



d'^r(|<5a>^|^ - IcJ*^!^) = 1, 



(68) 



This process occurs mostly at zero temperature, which 
is a pure quantum effect. The total damping rate is the 
sum of the two damping coefhcients: 7 = 75 + 7l- 

In this paper we are interested in homogeneous Bose 
gases, i. e. Ve-Kti^) — ^or homogeneous systems 
the condensate density remains the same throughout the 
space: $0 = y/TT^j while the excitations and the fluctua- 
tions satisfying Eq. ([50)1 and Eq. (|^ can be described 
by the plane waves 



(5$i(r) 
<5$2(r) 



1 



d'^ke''"' 



Uk 



Uk'(r) 
Wk' (r) 



d-'k', 



Vk' 



where Uk and Vk satisfy the Bogoliubov relations: 
_ie'ik)+gjnl)^/' + eik) 



ut = l 



2e(k) 



UkVk 



2e(k) 



(73) 



(74) 



(75) 
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FIG. 4: Landau damping rate per unit energy in two dimen- 
sions in the unit of Stt / 1 In 1 . The black solid line represents 
the result from HFB and the blue dashed line from the hy- 
drodynamical approach. The high-temperature asymptotic 
behavior and the low-temperature limit are also reported as 
the green dashed-dot line and the red dashed-dot-dot line, 
respectively. 



and e(k) is the Bogoliubov energy 



A. Quantum regime 

As mentioned in the last section, the dacay rate is 
mostly contributed by Bcliaev damping, which can be 
obtained by setting /° = / ° = O in Eq. The matrix 

element i?k'-k.k' reads 



B, 



k'-k,k' — 



[{uk{uk'Vk'-k + Vk'Uk'-k + Uk'Uk'-k) 



+Vk{uk'Vk'-k + Vk'Uk'-k + Vk'Vk'-k) 



(76) 



and other elements are zero. At zero temperature, only 
the momenta with long wavelength are involved in the 
Beliaev damping process, i. e. fc' ~ fc ~ |k' — k| <C mc. 
Therefore the long-wavelength approximation for the en- 
ergy e(k) and the wave functions Uk and Vk can be used: 



Uk 



e(k) ~ ck 



mc\i/2 1 f k 



/"mc\ 



fc3 
8^' 



1/2 



(77) 



1 / k 



8 V 2mc 



2 \2mc 

5/2 



1 / k 



2mc 



3/2 



Vk 



/"mc\ 
V 2fcy 



1/2 1 / fc 



1/2 



1 / k 



2 V 21110 J 8 V 2toc 



3/2 



1 / k 



8 v 2mc 



5/2 



(78) 



Substituting ^ and (|78D in Eq. 
result 



k,k' 



3 |k||k'||k'-k| 



V 4V2 {mc 



76l). we obain the 



(79) 



|3/2 



Inserting Eq. ([75]) to Eq. ([7^ and summerizing all the 
momenta k and k', one obtain the same result (|19p as 
that using hydrodynamical approach. Therefore we re- 
produces the results for 3-D and 2-D decay as Eq. pp)) 
and Eq. 



B. Thermal Regime 

In the theraml regime where luq ~ ksT, a long- 
wavelength Goldstone mode with the eigenfrequency 
(cjo — ck) describes the behavior of the condensate in 
the thermal clouds. In this limit, the u and v functions 
can be expanded as 



Uk 



Vk 



1/2 



1/2 



2e(fc)j 



Y 

2 \2mc^ J 
2e{k) J 2 V2mc2 J 



(80) 



Using this expansion,the long-wavelength behavior for 
the nonzero elements of the matrix A can be expressed 
as 



A 



k'. k' 



V \2mc'^ 



1/2 



Uk'Vk' 



'g '2-ul,vl, 



(81) 



with the angle O between the vectors k' and k, and the 
group velocity of the excitation defined as Vg — dek/dk. 

In three dimensions, the damping rate can be obtained 
by inserting the nonzero coefficients ([5T|) into ([7T|l and 
integrating out the angle 9 as foUows: 



.4=3 



(82) 



where r = ksT/mc^ as dimensionless temperature, 03 
is the three-dimensional scattering length, and F{t) is 
defined in the foUowing: 



(83) 



with the definition u = VT + z^. This expression was 
first found by Pitaevskii and Stringari [T^ . 

For low twmperature fesT <SC mc^, i. e. r <C 1, the 
function F takes its limit F « Stt^^^t^ /5 and one finds 
the Hohenberg and Martin's result (I^S)) . As mentioned in 
the last section, the hydrodynamical approach and HFB 
give the same limit at low temperature. However, at high 
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temperature, the hydrodynamic approach fails. For tem- 
perature r ^ 1., i. e. k^T 3> mc^, the function F takes 
the asymptotic Umit F Stt^^^t/A, and the damping 
rate approaches the Szepfalusy and Kondor result (^7)1 . 
Therefore the HFB gives correct asymptotic value at high 
temperature. 

In Fig. [3] the famous result for the three-dimensional 
damping rate using HFB method obtained by Strin gar i 
and Pitaevskii , and then recovered by Giorgini |13| 
is shown as the sohd Hne. We can see that the three- 
dimensional damping rate leaves the low-temperature 
limit very soon, while it approaches the high-temperature 
linear law very slowly. 

In two dimensions, HFB offers a good approximation 
for all regime of temperature. After inserting the matrix 
elements ([5T|) into [71] and then integrating out the angle 
6, one obtains the damping rate: 



e{k) e{k) 



2m52G(T), 



(84) 



where 



x/2 r , sech2(^) 



G(r) =^ / dz 



2(z^-l) 



(85) 



1 1 

2^ ^ 2^ 



In the low temperature limit r ^ 1, the function G takes 
the limit: G(r) ^/Sttt^ /16, the damping rate goes 
to the result As in three dimensions HFB gives 

the same result as the hydrodynamical approach at low 
temperature in two dimensions. 

For high temperature fc^T 3> mc^, the function G 
takes its asymptotic value: G(r) 0.013t, and the 
damping rate is given by 



7 



d=2 



0.026 



mkBT 



e{k) no 
Therefore the damping rate itself reads 



.,d=2 



7 



0.013, 



Stt 



kksT 



ln(l/na2)| mc 



(86) 



(87) 



In Fig. [4l the two-dimensional damping rate per unit en- 
ergy using HFB method is also shown as a function of r 
(solid line). We can see that the two-dimensional damp- 
ing rate approaches the high temperature linear law much 
sooner than the three-dimensional case. That means, 
the two dimensional systems go to the high-temperature 
asymptotic value at lower value of temperature compared 
to that in three dimensions. This behavior has been 
found by Guilleumas and Pitaevskii studying a quasi two- 
dimensional System (s^ . In this figure one can also see 
that the hydrodynamical approach can only give good 
results in the regime where t < 0.5. 



V. CONCLUSION 

In this work, we have compared the hydrodynamical 
approach and the HFB approach to calculate the damp- 
ing rate in 2D and 3D bose gas. The hydrodynamical 
approach is a powerful tool due to the fact that one can 
use the Green's-function technique based on the effective 
action ^ . This works very well at zero and low tempera- 
tures. However, this method truncates the rapid oscillat- 
ing fields, which is not the case for the high temperature 
regime, therefore it overestimates the Landau damping. 
On the other hand, the HFB approximation can explain 
either low temperature or high temperature damping. It 
seems that the mean-field approach (HFB) is a better 
method for Bose gases. The HFB approach based on the 
mean field method factorizes the quartic terms and there- 
fore Shohno Ansatz has to be used to avoid anomalous 
behavior. In the absence of the Shohno Ansatz, there 
would exist a gap in the low excitation spectrum. There- 
fore the mean field approach cannot guarantee a zero 
energy gap. From a physical point of view, the existence 
of a gapless excitation is a general rule for Bose systems. 
In order to avoid errors in the higher order calculations, 
the mean field approach should be very carefuUy used. 
Therefore we can see the benefit of the hydrodynamical 
approach for the low temperature regime. The low en- 
ergy excitation is always gapless using hydrodynamical 
approach. Another benefit of using the hydrodynami- 
cal approach, as indicated by Popov [lB|, is that it can 
avoid the strange divergence at high and low momenta, 
so-called ultraviolet and infrared catastrophe, which can 
be caused by the perturbation theory based on the mean- 
field approach. 

We have found for the first time that the Beliaev damp- 
ing rate is proportional to k^ at zero temperature and the 
Landau damping rate for the 2D bose gas is proportional 
to for low temperature and to T for high temperature. 
The behavior of the 2D damping is also totally different 
from the 3D damping. While the 3D Landau damping 
approaches the linear regime very slowly with increasing 
temperature, the 2D damping become linear very fast. 
The linear regime symbolizes the classical high tempera- 
ture behavior, therefore the two dimensional systems go 
to the high-temperature asymptotic value at lower value 
of temperature compared to that in three dimensional 
System. This behavior was also found in Ref. [32l | with 
numerical calculation for a quasi-2D system. 
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